Spin heat accumulation and its relaxation in spin valves 
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We study the concept of spin heat accumulation in excited spin valves, more precisely the effec- 
tive electron temperature that may become spin dependent, both in linear response and far from 
equilibrium. A temperature or voltage gradient create non-equilibrium energy distributions of the 
two spin ensembles in the normal metal spacer, which approach Fermi-Dirac functions through en- 
ergy relaxation mediated by electron-electron and electron-phonon coupling. Both mechanisms also 
exchange energy between the spin subsystems. This inter-spin energy exchange may strongly affect 
thermoelectric properties spin valves, leading, e.g., to violations of the Wiedemann- Franz law. 



The electric conductance through ferromagnet| normal 
metal|ferromagnet spin valves is a function of the mag- 
netic configuration.!. It reflects the spin accumulation, 
i.e., the spin (index a) dependent chemical potential \x a 
of the normal-metal island. The latter parameterizes 
the spin dependence of the energy distribution functions 
fa(E), whose description also requires spin-dependent 
temperatures T a ££ As shown below, these should in gen- 
eral be interpreted as effective parameters. 

In this Rapid Communication we describe the pro- 
cesses affecting the T a and through them the thermoelec- 
tric response in spin valves, which we find to be a sensitive 
probe for the non-equilibrium state in the non-magnetic 
spacer. Whereas the spin accumulation relaxes only by 
scattering processes that break spin rotation invariance 
such as spin-orbit interaction and magnetic disorder, the 
spin heat accumulation T s = T-j- — Tj, is sensitive also to 
electron-phonon (e-ph) and electron-electron (e-e) inter- 
actions. Spin-flip scattering in Al, Ag, Cu, or carbon 
is weak and hardly temperature dependent; the typical 
spin-flip scattering time r s f is of the order 100 ps, 2 which 
can be much longer than the dwell times in magnetoelec- 
tronic structures. The inter-spin energy exchange rate 
due to inelastic effects is strongly temperature depen- 
dent and above cryogenic temperatures typically dom- 
inates the direct spin-flip scattering in dissipating the 
spin heat accumulation. The spin heat accumulation in 
normal metal spacers should not be confused with the 
spin (wave) temperature of ferromagnets.— 

In a spin valve (Fig. [TJ , a nonmagnetic island is cou- 
pled to two ferromagnetic reservoirs with parallel (P) 
or antiparallel (AP) magnetic alignments. The chemi- 
cal potential of the left (right) reservoir is (J-l(r) and the 
temperature is Tt,(r). The conductances Gi/jj ff and See- 
beck coefficients St, / Ba of the contacts between the island 
and the reservoirs depend on spin a G {til}- Biasing 
the spin valve with either a voltage AV^ = (fiR — 
or a temperature difference AT = T B — T^ gives rise 
to a spin-dependent energy distribution function fa-(E) 
of the electrons on the island. As shown below, in the 
linear response regime this can be described exactly by 
spin-dependent chemical potentials and temperatures, 
such that f a {E) = f (E;n a ,T„), where fo(E;fX,T) = 




FIG. 1: (Color online): Schematic spin- valve biased with a 
voltage and/or temperature difference. Spin-flip and inelas- 
tic electron-electron and electron-phonon scattering in the 
normal metal spacer lead to inter-spin energy exchange and 
change the thermoelectric characteristics. / and Q stand for 
the charge and heat currents flowing into the island. T p h is 
the temperature of the phonon bath. 

{exp[(£ - (j,)/(k B T)] + is the Fermi-Dirac distribu- 
tion function. \i a and T a are determined by conservation 
of charge, spin and energy (see Eqs. ([2]))- The response 
matrix of the spin valve 
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relates the charge and heat currents / and Q to the biases 
AV and AT, respectively. Below, we derive expressions 
for the heat conductance K and thermopower S, in the 
presence of inter-spin energy exchange and for different 
magnetic configurations. 

The steady state potentials and temperatures can be 
determined from Kirchhoff 's laws for charge and energy 
for each spin. 4 For small eAV/k B , AT < T t , Tj., 

h,* + G si - M-a)/e = (2) 

i=L,R 

<2i, ff + K^(T a - T_ a ) + if - ph (T a - T ph ) = 0. 

i=L,R 

Here = G icr (^ CT -^)/e+G icr 5j cr (T cr -T l ) is the charge 
current for spin a through contact i, Qi tCr — Ci)Gi a T(T a — 
Ti) + Gia-Sia-T(^ a — fii)/e is the corresponding heat cur- 
rent, Gi a and Si a are the associated charge conductances 
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and Seebeck coefficients, and Cq = it 2 k 2 B / (Ze 2 ) is the 
Lorenz number. Spin decay is described by the (intcr- 
)spin conductance G s f = e 2 vpQ,/T s f for an island with 
volume fi, density of states at the Fermi level vp and 
spin- flip relaxation time r s f. The term A e _ p h describes 
the interaction with the phonons at temperature T p h. 
Inter-spin energy exchange is governed by the spin heat 
conductance = CoG s [T + K^_ e , where the first term 
originates from the spin-flip scattering and the second 
is due to e-e interactions. We are allowed to discard the 
spatial dependence of the distribution functions when the 
diffusion time Tp, = L 2 /D in the island with length L and 
diffusion constant D is shorter than both r s f and the spin 
thcrmalization time r st = Coe 2 vpTQ,/ '(A c _ p h + 2K<^). 

The in general lengthy solutions of Eqs. ([2]) are con- 
siderably simplified for left-right symmetric conductances 
and Seebeck coefficients, parameterized by Go = Gf+G^, 
P = (G t - GJ/Go, S = {G t S t + G^/Go and 
P' = (G t Sf - G^S^/iGoSo) for both junctions. In the 
antiparallel case the signs of P and P' in one of the junc- 
tions are inverted. In the parallel configuration the heat 
conductance becomes 



Kp = CqGpT 



2A c , ph r(l - P 2 7 ) 
1 _ P 2 7 + K c _ ph /(£ G T) 

and in the antiparallel configuration it is 

2A c _ ph r 



A AP = C G P T(1 - P 2 7 ) 



l + # e _ ph /(£ GoT)- 



(3) 
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The factor r = (T ph - T L )/(T R - T L ) ~ 1/2 parame- 
terizes the phonon temperature on the island: If the 
phonons are poorly coupled to the substrate, as for ex- 
ample in perpendicular spin valves or in suspended struc- 
tures, T ph = (Tf + T i )/2. For the P configuration this 
yields r = 0, whereas for the AP configuration we get 
r = -A ti P/[2(A c _ ph + K n + C G T)]. In the op- 
posite limit r = ±1/2, viz. T p h is fixed to the bath 
temperature of the left or right reservoir. The coeffi- 
cient 7 = [1 + (AVph + 2K U )/(£ G T)}- 1 describes 
inter-spin energy exchange. Factoring out the temper- 
ature dependence of A e - P h oc T 4 and K^_ e cx T v+1 
(see the discussion below) yields 7 = [1 + (T/T c h, p h) 3 + 
(T/Tch.o-c)^ + 2G s f/Go] _1 , where the characteristic 
temperatures are T ch ,o- P h = [(£oG T 4 )/AVph] 1,/3 , 
T ch , e -e = [{C Q G Q T' J + 1 )/{2Kli c )] 1 /' J for elcctron-phonon 
and electron-electron couplings, respectively. The expo- 
nent v depends on the dimensionality (nd) of the sample. 
We are here mainly interested in 3d samples (y = 3/2) in 
which all sample dimensions exceed the thermal coher- 
ence length £ T = v / hD/{2nk B T). 

In the parallel configuration the thermopower satisfies 
Sp = Sq and in the antiparallel one£ 

S AP 1-PP' + 2G sf /G + 7 P (P - P' - 2P'G si /G ) 



Sp 



1 - P 2 + 2G sf /G 



(5) 

The temperature dependence of A and S is plot- 
ted in Fig. H for A c _ ph > A c _ c , C G sf T. For T < 
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FIG. 2: (Color online): Temperature dependence of the heat 
conductance K (solid lines, left axis) and thermopower S 
(dashed line, right axis) of a structurally left-right symmetric 
spin valve with P — 0.9, P 1 = 0.5 and when the electron- 
phonon relaxation dominates the inter-spin energy exchange. 
The lines are plots of Eqs. I©-© and the symbols have been 
calculated from the full nonequilibrium distribution function, 
Eqs. (|10|) and pip . The results have been calculated for P 
configuration with r = (circles) and r = 1/2 (squares) and 
AP configuration with r — (stars) and r = 1/2 (triangles). 



min(T c h, c -c, T c h, c -ph) = T c h, the device operates as a 
spin heat valve in which the heat current can be con- 
trolled by the magnetization configuration. Contrary 
to the charge conductance, however, the magnetoheat 
conductance (Kp — Aap)/Ap vanishes for T ^> T c h or 
7^0. Thus the presence of inelastic scattering leads 
to a violation of the Wiedemann- Franz law A = CqGT 
for T > T c h. The magnetothermopower (5p — Sap)/Sp 
persists provided P ^ P'.— The measured heat conduc- 
tance and thermopower as a function of temperature and 
magnetic configuration may hence yield unprecedented 
information on the energy relaxation in normal metals. 

We now address the characteristic temperatures 
Tch,c-ph and T c h, c -c- The former can be obtained directly 
from the Debye form for the heat conductance between 
electrons and acoustic phonons^ A c _ p h = |Xf2T 4 , valid 
for T <C Tbcbyc- Here E is the e-ph coupling constant® 
and the factor 1/2 takes into account spin degeneracy. 
The characteristic temperature for electron-phonon cou- 
pling thus reads 



T, 



ch,' 



- P h 



I 157i£fi J 



1/3 



G h 



1/3 



(6) 



For T > Toobyo, the electron - acoustic phonon scat- 
tering and thereby inter-spin energy exchange saturates. 
Optical phonons start to contribute in this temperature 
regime, but are disregarded here. 

The e-e scattering collision integrals with spin- 
dependent distribution functions contain three terms 



A-c,o-(e) 



T(7(7 

J (a) 



•(e) 



presented by the diagrams in Fig. |3J 

Processes (b) and (c) induce inter-spin energy ex- 
change, which can be described in terms of a heat current 



3 



(a) a 



(b) a 



(c) 3 



e + Tiu)\^ fojj I e' — Hu € + huj\^ f luJ j e — Til: e + Tilj\^ ^ j e —TkjJ 
»N/v\/» ^\yx/\* • r xf\/\s 



FIG. 3: Electron-electron scattering vertices, (a) Equal-spin 
scattering, which equilibrates the electrons but does not ther- 
malize the spins, (b) Spin conserving scattering and (c) spin 
exchange scattering, which do thermalize the spins. 



flowing between two spin ensembles, 



dee(lj l 
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The direct spin current due to e-e interaction vanishes 
in the absence of spin-orbit scattering, J de(lj^ 



V) 



I(X) — 0. In 3d, to lowest order in spin particle and heat 
accumulation, — /i^) / (/i^ + fi^) <C 1 and T S /(T^ + 



n) « i, 

where 

K n _ 



we arrive at Q, 



n 



105C(7/2)4 /2 ^ 5/2 
32[2ttEt{1 + F)f/ 2 h' 

FC [tt 2 C(3/2) + f|C(7/2)]4 /2 T 5 / 2 



(8a) 



(8b) 



K (c) = (F + 2) 2[2ttE t (F + \)f/ 2 h 

Here C = -1 + (F + l) 3 / 2 , F > -1 is the spin triplet 
Fermi liquid parameter (F = — 1 corresponds to the 
Stoner instability), Et — hD/Q 2 / 3 is the Thouless energy 
proportional to the inverse time it takes to diffuse over a 
length SI 1 / 3 and £(x) is the Zeta function. Summing the 
two contributions from Eqs. ([5]) yields the characteristic 
temperature 

r ch , c -c = 8 7 r^V /3 T^(^ + l)x 0) 



tt(F + 2) 



48 j FCtt 2 C(3/2) + 105[6 + F(3 + G)]C(7/2) 
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In Id and 2d structures the spin-flip contribution (c) 
has an infrared divergenc o 9 ' 10 that needs to be regular- 
ized. As a result, the inter-spin energy exchange due to 
e-e scattering becomes stronger and the corresponding 
7ch.c-c l° wer - This is especially relevant at low temper- 
atures and small structures since £t may exceed 100 nm 
at T ss 1 K. We intend to analyze the resulting inter-spin 
energy exchange in reduced dimensions in the future. 

In order to assess the relevance of our results for re- 
alistic samples we consider a disordered island of a spin 
valve coupled to the reservoirs via tunnel contacts. For 
example, with F = —0.3 we get 

2/3 
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Making the sample smaller and conductance larger in- 
creases both characteristic temperatures, but the increase 
for T c h. -ph is slower. For ft = 0.001 (/im) 3 and Go = 1 
S we get T ch , e -ph = 22 K whereas T cr , _ c = 400 K. We 
may therefore conclude that in spin valves with metallic 
contacts and 3d spacers the inter-spin energy exchange 
due to e-e interaction can be neglected. The spin ther- 
malization rate with F = —0.3 is 



1 

1 

25ns 



1 



20ns 
T 

Ik 



T y /2 ^0.001 m 2 /s\ 3/2 
IK, 



D 



10 9 Wm~ 3 K 



-5 



10 4 



J" 1 m- 3 



Up 



The first term comes from e-e scattering and the sec- 
ond from e-ph scattering. This rate exceeds the spin-flip 
scattering rate ~ 10 GHz at temperatures above ~ 10 K. 

Above we assume that the electron energy distribution 
function is well represented by Fermi-Dirac distributions 
with spin-dependent chemical potentials and tempera- 
tures. This is not true in general, since / ff (e) has the 
nonequilibrium form^^i 



GlctIl + G Ra f R + is F e 2 m con [f a , f_ a ] 



G 



La 



Gb.o 



(10) 



where fi,m = /o(e; ^LIRi T) are the distribution func- 
tions for the reservoirs and / co ii describes all inelastic 
scattering events. The charge (n — 0) and heat (n = 1) 
currents through contact i then become 



Gi, 



(e)(U(e)-fi(e)). (11) 



Thermoelectric effects can be included by adding a 
weak energy dependence to the conductances, Gi a {t) rs 
Gia [1 + c ia (e - Mi)]; and expanding to linear order in Ci, a . 
Identifying — e£oc ia T, we recover Eqs. (01 and §5§ 
in the regime eAV/k B ,AT < T L ,T R sa T even in the 
absence of collisions (i.e., 7 = 1). For Cj CT = and to 
linear order in the applied bias, the nonequilibrium dis- 
tribution pop is identical to the quasiequilibrium one. 
Under these conditions, the collision integrals can be cal- 
culated by replacing the full distribution functions by the 
quasiequilibrium ones. Numerical solutions of the kinetic 
equations (see Fig. [2} indicate that in linear response col- 
lisions and finite Ci CT 's do not change this conclusion. 

Beyond linear response spin-dependent temperatures 
can strictly speaking be invoked only in the presence 
of strong inelastic scattering such that Tf s» Tj.. Nev- 
ertheless we can define effective electron temperatures 
that satisfy the standard relation with the thermal en- 
ergy density in the Sommerfeld expansion^ 



T 



-TT-J r ^ ~l + e(e- ^)}ede. (12) 



Proceeding with Fermi-Dirac distributions with effective 
spin-dependent temperatures and chemical potentials, \i a 
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FIG. 4: (Color online): Spin-dependent effective tempera- 
ture vs. voltage in an asymmetric spin valve with P = 0.9, 
P' = 0.5 and Gr = O.lGz,. The lines are calculated from 
Eqs. @ and ((15)) and the symbols from Eq. (TT!^ for nu- 
merical solutions of the kinetic equations. The upper curves 
are for majority, the lower for minority spins, and different 
strengths of e-e scattering with F = 0: no scattering (solid 
line and circles), weak scattering with Et = 0.05fcsT and 
Gl — 100e 2 //i (dashed line and squares), and strong scatter- 
ing with E T = 0.001fc s T and G L = W0e 2 /h (dash-dotted line 
and stars) . Here T denotes the temperature of the reservoirs. 
Left inset: behavior at low bias with thermoelectric effects 
cl = IOcr = 0.02/(fcT). Right inset: distribution function at 
eAV = IOEt with different strengths of e-e scattering. 

and T a can be obtained from Eqs. @ by replacing the 
expression for the charge and heat currents through con- 
tact i with their nonlinear counterparts, 

I a = ^ { Mct - m + ^ [£ oe 2 (T 2 - 2?) - ( Pa - ^f] } 

g,, CT = G i<a [£ (T 2 - T?)/2 - (,4 - ^)/{2e 2 )] (13) 
+ GW(/V - w ) [£ (T 2 + Tf)/2 - - ^ 2 )/(6e 2 )] . 

These equations are obtained by a direct integration of 
Eq. (fTTj) using Fermi-Dirac functions fi(e) and /o-(e). We 
also have to replace the linear-response forms of the spin 



mixing terms in Eqs. ([2]) by their forms far from equilib- 
rium. For example, for e-e scattering with F — we use 
Q aS = 15C(7/2)4 /2 (rJ /2 - Tl^)/[16n(2%E T f/ 2 ]. 

In the absence of collisions and for weak thermoelec- 
tric effects it can be proven by direct integration that 
the effective temperatures defined by Eq. (fT2|) agree with 
those which follow from heat conservation. In Fig. |4] we 
present a complete numerical solution of the kinetic equa- 
tions along with the results from the quasiequilibrium 
heat balance equations from which we conclude that the 
two approaches for calculating T a agree also in the pres- 
ence of inter-spin energy exchange. 

Spin heat accumulation cannot be directly measured 
by two-terminal transport experiments in linear sys- 
tems. In order to prove the presence of a sizable T s far 
from equilibrium it should be probed by spin-selective 
thermometry, such as a generalization of the tunnel- 
spectroscopy in Ref. 1 1 lL by measuring the shot noise of 
the spin valve, or through electron spin resonance. 

In conclusion, we have shown that inter-spin energy 
exchange in a spin valve affects the temperature and 
magnetic configuration dependence of its thermoelectric 
properties. The different thermalization mechanisms can 
be quantified by characteristic temperatures, Eqs. ([6]) 
and ([9]), above which interaction effects become impor- 
tant. We introduce the concept of spin heat accumu- 
lation via the spin-dependent effective electron tempera- 
tures T a in Fermi-Dirac distribution functions, which can 
be used to describe transport properties beyond the lin- 
ear response regime. We demarcate the regime in which 
spin valves can be employed to control heat currents. 
Other types of operations can be envisaged as well, such 
as spin-selective cooling of the electrons (see the left inset 
of Fig. HJ). 
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